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Complexation of a long flexible polyelectrolyte (PE) molecule with oppositely charged spherical 
particles such as colloids, micelles, or globular proteins in a salty water solution is studied. PE binds 
spheres winding around them, while spheres repel each other and form almost periodic necklace. If 
the total charge of PE in the solution is larger than the total charge of spheres, repulsive correlations 
of PE turns on a sphere lead to inversion of the net charge of each sphere. In the opposite case 
when the total charge of spheres is larger, we predict another correlation effect: spheres bind to 
the PE in such a great number that they invert the charge of the PE. The inverted charge by 
absolute value can be larger than the bare charge of PE even when screening by monovalent salt is 
weak. At larger concentrations of monovalent salt, the inverted charge can reach giant proportions. 
Near the isoelectric point where total charges of spheres and PE are equal, necklaces condense into 
macroscopic bundles. Our theory is in qualitative agreement with recent experiments on micelles-PE 
systems. 

PACS numbers: 87.14Gg, 87.16.Dg, 87.15.Tt 



I. INTRODUCTION 



Electrostatic interactions play an important role in 
aqueous solutions of biological and synthetic polyelec- 
trolytes (PE). The complexation of a long flexible poly- 
electrolyte with appositely charged pSpherical particles 
such as micelles, El globular proteinsB or colloidsB is a 
generic electrostatic problem of the polymer physics. A 
long PE binds oppositely charged spheres winding around 
each of them (Fig. 1). Without losing the generality, we 
assume that the PE is negative and spheres are positive. 

If the charge of a sphere is not completely compensated 
by the winding PE, the net charge of the sphere is still 
positive, the neighbouring spheres repel each other and 
form on the PE an almost periodic necklace (Fig. 1). The 
same picture is true when the winding PE inverts the net 
charge of each sphere making it negative. We call this 
nontrivial phenomenon a sphere charge inversion (SCI). 
SCI is known to happen in the most famous biological 
example of PE-spheres complexation. In the chromatin, 
the negative double-helix DNA molecule winds around 
a positive histone octamer to form a complex known as 
the nucleosome bead. Nucleosome beads are connected 
by DNA linkers in the so-called beads-on-a-string struc- 
ture. When linkers are cut enzymatically each nucleo- 
some bead is found to have a negative net charge. 

The counterintuitive phenomenon of SCI has attracted 
a lot of attention of theorists. However, all theoretical 
and numerical studies of SCI, have been done fbr_Jjhe 
complexation a single sphere with a PE molecule □113. 

In this paper, we propose the first theory of the SCI in 
the necklacelike complex of the PE with many spheres. 
Our theory accounts for the interaction between different 
spheres. Wepajrgue that in this case, as in the case of a 
single spheret2l, SCI happens due to repulsive correlations 
of different PE turns on the surface of spheres. 



For many spheres, however, not only the net charge of 
a sphere should be found, but simultaneously the number 
of spheres attached to PE molecules is to be calculated. 
Therefore, the second and even more challenging problem 
is to determine the sign of the whole complex of PE with 
many spheres. Is it positive or negative at given number 
concentrations of PE, n^, and spheres, n^, in solution? 

The standard Debye-Hiickel and Poisson-Boltzmann 
theories of screening of PE by monovalent counterions 
leave the net charge of PE always negative. These the- 
ories, however, do not work for screening by strongly 
charged spheres which, as we mentioned above, form 
a correlated sequence, reminding a necklace or a one- 
dimensional Wigner crystal. One can call it a Wigner 
liquid, because the long range order in many practical 
situations is destroyed. 
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FIG. 1. The beads-on-a-string complex of a negative PE 
molecule and many positive spheres. On the surface of 
each sphere, due to the Coulomb repulsion, neighbouring PE 
turns lie parallel to each other. Locally, they resemble an 
one- dimensional Wigner crystal with the lattice constant A. 
At a larger scale, charged spheres repel each other and form 
another one-dimensional Wigner crystal along the PE with 
lattice constant x + 2R. A Wigner-Seitz cell of this crystal is 
shown by the thick arrows. 

Wigner-crystal-like correlations between multivalent 
count eriorp,ax^inown to lead to charge inversion of rigid 
macroionsEjEj'Q. This happens because when a multi- 
valent ion approaches an already neutralized macroion. 
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it repels other counter ions, creates for itself a correla- 
tion hole or an image of opposite charge which attracts 
it to the surface. In the necklace shown above, the PE 
segment wound around each sphere interacts exclusively 
with this sphere and plays the role of the correlation hole 
or a Wigner-Seitz cell. Therefore, again, Wigner-crystal- 
like correlations come into play and lead to an additional 
attraction of the spheres to the PE. Indeed, when a new 
sphere approaches a neutralized necklace, it pushes other 
spheres away, unwinds a segment of PE from them and 
winds this segment around itself. This segment is the 
sphere's correlation hole or, in other words, its image 
in the PE. We are dealing with the correlation physics 
because the image appears only in response to the new 
sphere. We show below that, at large rig and small n^, 
this correlation attraction leads to PE charge inversion 
(PECI). PECI was observed in a micelle-PE systemEJ. 

Following Ref. pT[-p!3| for a quantitative characteristic 
of charge inversion we introduce the charge inversion ra- 
tio of the PE V = —Q'^/Q, where Q = Lr] is the negative 
bare charge of PE {L and 77 are, respectively, the contour 
length and the linear charge density of a PE molecule) 
and Q* is its positive net charge together with all ad- 
sorbed spheres. Optimization of the free energy of a 
complex with respect of the number of bound spheres 
per PE molecule. A/", shows that, even for a large Debye- 
Hiickel screening radius of the solution, the optimal 
N = No is so large that 

Here R and q is the radius and charge of a sphere and a is 
a dimensionless logarithmic function of qR/rjr'j) (see Sec. 
III). We assume everywhere in this paper that q/r]R ^ 1, 
so that more than one turn of PE winds around the 
sphere to neutralize it. 

PECI also grows with stronger screening (smaller td)- 
For in the range A <^ ro <^ R{q/ Rr]Y^'^ ^ we show 
that 

Here /3 is a dimensionless function of q/Rrj and R/vd 
(see Sec. V). A PECI given by Eq. (0) and Eq. (|) can 
be called giant. 

At a large sphere concentration n^, when the PE num- 
ber concentration, n^, grows and reaches ng/TVo, the pool 
of free spheres gets exhausted and each PE molecule can 
not get the optimal number Nq of them any more. Then 
PECI becomes weaker and disappears linearly at the iso- 
electric point Upi = qns/\Q\^ where the total charge of 
all the spheres compensates the total charge of all PE 
molecules. When the concentration Up continues to grow 
beyond Upi practically all the spheres remain bound to 
PE and the net charge Q* is negative and grows by ab- 
solute value. Variation of (5*/(5 with rip is shown on Fig. 
2 by the solid line. 





1 


Q*/Q 




r 




h 














q*/q 


PECI / 




SCI 















p 



FIG. 2. Schematic plot of q jq and Q* /Q as functions of 
the PE concentration, rip, at a fixed and large sphere concen- 
tration ris and at a not very large td- The shaded stripe cor- 
responds to the region around the isoelectric point rip — ripi 
where necklaces condense into macroscopic bundles. 

Simultaneously with these variations of Q*/Q, the net 
charge q* of a sphere changes, too. At rip < rig/No it is 
positive and close to q. At > n^/A/'o, the net charge q^ 
starts to decrease linearly with rip — ripi. At the isoelectric 
point rip = ripi, the charge crosses zero and simulta- 
neously, the linker length x vanishes. At rip > n^^, the 
charge q^ becomes negative and SCI appears. We show 
that the charge inversion ratio of a sphere, S = —q*/q^ 
grows with rip — ripi until it reaches the value xiorrespond- 
ing to a single sphere bound to infinite PFHj, which is 
roughly equal to the inverse number of turns necessary 
for PE to neutralize a sphere. The behavior of q* / q as 
function of rip is shown by dashed line in Fig. 2. It is 
clear from Fig. 2 that SCI happens at rip > ripi and 
PECI at rip < ripi. 

Thus, we arrive at the conclusion that both at rip > ripi 
and np < n^^, a beads-on-a-string structure can sponta- 
neously self-assemble from a PE and oppositely charged 
spheres without any non- Coulomb forces. The latter 
structure resembles the lOnm fiber structure of the chro- 
matin. 

Experimental observation of SCI is possible when 
spheres with winding PE are cut out from the complex. 
Then their charge can be measured by electrophoresis. 

Consequences of PECI are more pronounced. PECI 
leads to reentrant condensation of necklaces into macro- 
scopic bundles. Indeed, near the isoelectric point rip = 
ripi each complex is almost neutral and short range dX^ 
tractive forces between Wigner-crystal-like complexestSi 
lead to their condensation and coaservation. Away from 
the isoelectric point each necklace complex is charged 
and their long range repulsive interactions prevent their 
condensation. One can watch how condensation begins 
and ends changing one of concentrations. For example, 
if we keep the spheres concentration ris large and fixed 
and start from rip ^ ripi, the complexes are negative 
and repel each other. Then with decreasing rip the con- 
densation happens in the vicinity of the isoelectric point 
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(the shaded region in Fig. 2). If we continue decreas- 
ing Up, FECI begins and the complexes become positive. 
When their positive charge, Q*, becomes large enough, 
the coaservate dissolves. An important prediction of such 
theoryllZl is that the electrophoretic mobility changes sign 
in the coaservation range. We estimate the width of the 
range of rip around ripi where coaservation occurs. This 
width increases with decreasing td and dit vd <^ A: 

ripi \ q J td ' 

The narrow range of coaservation followed by ijesolu- 
bilization was observed in the micelles-FE systemEl as a 
function of the charge of micelles. The electrophoretic 
mobility of complexes was indeed found to change sign 
within the interval of the micelle charge in which coaser- 
vation happens. The width of the coaservation region 
was also observed to increase with decreasing rjj in qual- 
itative agreement with Eq. (||). 

To illustrate the physical picture discussed above we 
carry out Monte-Carlo simulation of the complexation of 
a negative FE with two positively charged spheres. The 
system is in a salt free solution. The simulated spheres 
have charge 70e uniformly distributed over their surface 
and radius 3.51b, where Ib = 7. 2 A is the Bjerrum length 
at the room temperature. The FE is modeled as a chain 
of free jointed hard spherical beads with radius 0.21b and 
charge — e. The bond length is kept fixed and equal Ib- 
The Monte-Carlo algorithm is described in our previous 
paper (Ref. |l0|). 

The snapshots of three such complexes are shown in 
Fig. ^. In the first simulation, the FE molecule has 70 
monomers. This complex illustrates the regime where 
the spheres are in abundance {up < Us/Nq). In the 
second simulation, the FE molecule has 140 monomers 
so that the complex is neutral and illustrates the FE- 
spheres complexes near the isoelectric point Up = Upi. In 
the last simulation, the FE molecule has 210 monomers. 
This complex illustrates the regime where there are not 
enough spheres to neutralize the FE {up > Upi). 




FIG. 3. Snapshots of three complexes of a negative PE with 
two positively charged spheres. The numbers of monomers of 
the PE in the cases (a), (b) and (c) are 70, 140 and 210 re- 
spectively. All the spheres have charge 70e. The total charges 
of the complexes are 70e, and — 70e respectively. 

FECI is clearly observed in the first simulation. One 
FE with charge -70e complexes with the two spheres with 
charge 70e each, making a giant 100% FECI. Around 
the isoelectric point, the distance between the surfaces of 
the two spheres is practically zero (less than the length 
of one FE bond Ib)- Many such FE-spheres complexes 
condense into a large bundle around the isoelectric point. 
Far beyond the isoelectric point, the FE-spheres complex 
is stretched again. SCI is observed with around 85 FE 
monomers bound to each sphere (~20% SCI). 

Although a perfect solenoid conformation of FE is not 
observed in Fig. 3, one can clearly see that FE segments 
of different turns stay away from each other and locally 
resemble a one dimensional Wigner crystal which helps 
to lower the energy of the system. Globally, thermally 
excited soft bending modes with characteristic length R 
melt the solenoid into a compromised "tennis ball" con- 
formationtj. The difference in energy between a "tennis 
ball" and a solenoid, however, is small compare to the in- 
teraction between the spheres and the FE. This explains 
the agreement between the observed physical features of 
the simulated finite temperature systems and those pre- 
dicted by our zero temperature theory. 

Here, we jalso would like to mention recent Monte- Carlo 
simulationdl3 of complexation of a FE molecule of given 
length with many oppositely charged spheres. Results of 
this work are in qualitative agreement with Fig. 2. How- 
ever, we cannot compare them with our theory quantita- 
tively because in these simulations the parameter q/Rrj 
is not large. 

It should also be noted that the behaviour of charge in- 
version for FE-spheres complexes described above is qual- 
itatively similar to that of lipid-DNA complexes studied 
in Refs. |l^,^ where one sees both kinds of charge inver- 
sion as one moves away from the isoelectric point. 

This paper is organized as follows: In Sec. II, we study 
the free energy of the system and derive equations for the 
equilibrium values of x and N. In sections HI and IV, 
these equations are solved for the weak screening case 
where the screening radius r^^ is larger than the necklace 
period x + 2R. In Sec. V, we discuss condensation and 
resolubilization near the isoelectric point. Section VI is 
devoted to the strong screening case, <C x + 2R. We 
show that in this case FECI is much stronger. In Sec. 
VII, we derive the charge inversion ratio for a stiff (rod- 
like) FE and compare it to the result for an intrinsically 
fiexible FE obtained in previous sections. It is shown 
that at weak screening, FECI is stronger for fiexible FE 
while at strong screening, FECI is stronger for rigid FE. 
This means that if one stretches the FE by external force, 
some spheres leave the FE in the weak screening case and 
condense on the FE in the strong screening case. Finally, 
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in the conclusion, we discuss several important assump- 
tions used in this work. 



II. OPTIMIZATION OF THE COMPLEX 
STRUCTURE. WEAK SCREENING CASE 

Let us start by writing down the free energy of the com- 
plex of a PE with length L and charge density r] wind- 
ing consequently around N oppositely charged spheres of 
charge q and radius R (see Fig. 1). First, we assume the 
complex is in a low salt solution so that the screening 
radius r/:) is larger than the distance between two neigh- 
bouring spheres x + 2R. We call this situation the weak 
screening case. Taking into account that the length of the 
PE segment that winds around each sphere is {L/N — x)^ 
we have 



F{N,x) = 



Q 



* 2 



DN{x + 2R) 



In M + Nf{x) 



(4) 



where 



q 



* 2 



2RD 



2(7*77, x^2R 
■ m 



D 2R 

+ (- + 2i^)^ln^-(x + 4i^)^ln^ 



-^(L/N -x)^\n- 
D a 



T] , X 

X-— m — . 

D 2a 



(5) 



Here D is the dielectric constant of water. At a length 
scale greater than its period x + 2R^ the complex is a 
uniform rod of length N{x + 2R) and charge density 
Q*/N{x^2R). The first term in Eq. (|) is the self-energy 
of this necklace (the macroscopic self-energy). The loga- 
rithmic divergence of this energy is cut off at small dis- 
tances hy X -\-2R and at large distances by the length 



J\f{x^2R) = mm{rD,N{x^2R)} 



(6) 



where N(x + 2R) is the rod length. In the second term 
of Eq. (Q), f{x) accounts for the total energy of one pe- 
riod of the necklace. It is calculated as the energy of a 
Wigner-Seitz cell consisting of a sphere with two PE tails 
of length x/2. The first terms in Eq. (^ accounts for the 
self-energy of the adsorbed sphere with net charge at 
the PE. The second term accounts for the interaction of 
the sphere with the tails, the third and fourth terms ac- 
count for the interaction between the tails. The fifth and 
sixth terms are, respectively, the self-energies of the PE 
wound around the macroion (which is screened at dis- 
tance A between turns) and of the two straight tails with 
length x/2. It should be noted that writing down the sec- 
ond of Eq. (^ as Nf{x) we have neglected the difference 
between the end spheres with those in the middle of the 
PE. This is justified for a reasonably large value of N. 
It should also be noted that we neglected the entropy of 
the PE monomers in the tails and at the spheres surface. 
This is justified because the charge of the sphere is large 



and Coulomb energy is much larger than the thermal en- 
ergy of PE. 

As we will see later, when Up is away from the isoelec- 
tric point 77pi, the linker length x is much larger than R. 
This helps to simplify Eqs. (||) and (||). Approximating 
A R^ /{L/N — x) and keeping only terms of the highest 
order in the large parameter x/R^ one can rewrite these 
equations as 



^2 

F{N,x) = —N\nAf^Nf(x) 

X 



(7) 



{S + xf 
2R 



2{S + x)ln- 
R 



/-r N, (L/N — x)a , X 



where we introduce the PE length needed to neutralize 
one sphere C = q/r] and 



S = jC-L/N = QyNr] 



(9) 



so that (7* = r]{5 + x). From now on, we also write the 
energy in units of rj^ / D (hence, the energy has dimen- 
sionality of length) . 

At a given the optimal distance x can be calculated 
by minimizing the free energy F{N^ x) with respect to x. 
This gives, to the leading terms. 



dx 



InA/"- 



S^x 
R 



In- 



R 



. L/N-x ^ 
In = . 



R 



(10) 



The physical meaning of each term in Eq. ( p^Oj ) is quite 
clear. When one brings a unit length of the PE from the 
sphere surface to their tails, thereby increasing x, the 
four terms of Eq. ( [lo| ) are, respectively, the lowering in 
the system's macroscopic energy (with increasing x), the 
potential energy cost due to the attraction of the PE to 
the sphere, the potential energy gained due to the repul- 
sion of two PE tails of each sphere and finally the cost in 
the correlation energy at the surface of the sphere. This 
last term - the correlation energy term - needs further 
clarification. If the PE turns around a sphere were ran- 
domly oriented, its self-energy per unit length would be 
\n{R/a). In reality, due to strong lateral repulsion be- 
tween different PE turns, they lie parallel to each other 
and locally resemble a one-dimension Wigner crystal. In 
this ordered state, the self-energy per unit length of the 
PE turn is screened at distance A instead of R. This gives 
the energy ln{A/a) per unit length of the PE. The low- 
ering in the self-energy of the PE segment wound around 
a sphere (with length {L/N — x)) in the ordered state 
as compared to the randomly oriented state is equal to 
{L/N - x)[ln{R/a) - ln{A/a)] = {L/N - x) ln{R/A) 
{L/N — x) ln{{L/N — x)/R) and is called the correlation 
energy. The fourth term of Eq. ([1q| ) is its derivative with 
respect of x. A more detail discussion of this correlation 
effect can be found in Ref. ffol 
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In principle, one can solve Eq. ( |Tq|) numerically for 
X as a function of N and other parameters of the sys- 
tem L/R and jCJR. After that, one can substitute x{N) 
back into Eq. (0) and find the optimal value Nq from the 
equation: 



dF{N,x{N)) 



dN 



(11) 



N=No 



where /i^ is the chemical potential of spheres in the bulk 
solution. 

If the PE concentration is small (n^ < Us/Nq) then 
No and x{No) define the configuration of the complex. 
However, in the case rip > Us/Nq there are less than 
A^o spheres for each PE. In this case, N = Us/rip (with 
an exponentially small correction) and x{N) defines the 
configuration of the complex. 

Let us now study asymptotic limits in which Eq. (p!o| ) 
can be solved analytically providing clear physical picture 
of our system. 



III. A SINGLE POLYELECTROLYTE MOLECULE 
IN CONCENTRATED SOLUTION OF SPHERES. 
WEAK SCREENING CASE 

In this section we consider the case where the PE con- 
centration is small, Up < Us/Nq^ so that the optimization 
of F{N, x{N)) with respect to N is needed to get the op- 
timal configuration of the complex. Here and everywhere 
in this paper we assume the bulk sphere concentration is 
high enough so that one can approximate /ig = /2R 
(the self-energy of a bare sphere) neglecting the entropic 
part of the chemical potential. Equation ( |ll| ) can be 
rewritten as 



dF _^ 
2R~ dN ~ X 
{8±xf 
2R 



1 



2L Nx' 
2L 



\njV- 
2Nx' 



N{8^x) 8^x 



- (2£ + X + Nx') In ^ + (x + Nx') In ^ (12) 

2R R 

where x' = dx/dN . 

To solve Eq. ([lo| ) for x, we assume that C ^ RhiM 
or, in other words, the screening length is smaller than 
an exponentially large length, <C xexp(£/i?). As we 
see below, in this case 5 ^ x and the last two terms in 
Eq. ([iq|) can be neglected. This gives 



= 5^^^{R\nM)^ 



/2 



(13) 



Substituting Eq. (13) into Eq. ( |12[ ) and keeping only the 
highest order terms one obtains the equation 

{5R\nMf''^{25 + ZL/N) - {L/Nf/2 = , (14) 

which has consistent solution only if 5 ^ L/N. In this 
case, 



L f L/N \ 
hJ/J 



1/3 



and the solution for A^o = L/{x + 2R) is 
L L 



No 



(53/4(i?lnAA)i/4 - £3/4(^in^)i/4 ' 
The corresponding charge inversion ratio is 

Q* nqs f c y/'' 



V 



RlnM) 



> 1 . 



(15) 



(16) 



(17) 



From Eq. (^ and (|T^) , it is easy to see that the relative 
order of all the lengths in the system is 

5 L/N = /:^/^(i?lnAr)^/^ 

> X = £i/2(i?lnAr)i/2 ^ ^1^^ ^ ^ (^g) 

This order is consistent with the assumptions we started 
with. 

As we saw above, the two last logarithmic terms in Eq. 
( p!2|) are negligible. Therefore, the main driving force be- 
hind PECI is the gain in the self-energy of a sphere when 
PE winds around it reducing its net charge. It is the dif- 
ference between the left hand side and the second term 
of Eq. (p!2|). In other words, the sum of the self-energies 

^ /2RD decreases when PE distributes itself over larger 
number of spheres. This correlation effect overcomes the 
macroscopic energy cost of overcharging the PE (the first 
term on the right hand side of Eq. (p!2|)). Therefore, 
PECI can be well obtained in the approximation where 
PE charge is smeared on the surface of spheresBH. 

Let us explain why we still call PECI calculated here 
a correlation effect. As we saw above the reason for this 
PECI is that each sphere is bound to several turns of a 
negatively charged PE. These turns can be considered as 
a correlation hole in the sense that this is the part of PE, 
which interacts almost exclusively with the given sphere 
(other spheres are at much larger distance x ^ R). The 
segments of PE wound around each sphere have the same 
length, L /TV — x ~ L/N ^ x. Therefore, similarly to 
Ref. |1Q|-P^ we are dealing with Wigner-crystal-like cor- 
relations and the wound segment can be considered as 
a Wigner-Seitz cell of the bare sphere. The gain in the 
sphere self-energy mentioned above is nothing but the 
usual binding energy per sphere of a Wigner crystal: the 
interaction of a sphere with its Wigner-Seitz cell. 

Note that because most of the PE length is wound 
around the spheres, the periodicity of positions of spheres 
covered by PE solenoids in the real space (see Fig. 1) is 
less important than in the case of a rigid PE (see Sec. 
VII) or other cases of chargeJrpLersion of rigid macroions 
by multivalent counterionsEJilj. Linkers between dif- 
ferent pairs of neighbouring spheres may differ in their 
length without a substantial change in the two major con- 
tributions to the free energy discussed above (the sphere 
self-energy gain and the macroscopic charging energy). 
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The thermal motion can even melt the Wigner crystal of 
spheres in the real space while the length of the wound 
segment remains unchanged. Therefore PECI is much 
more robust than the Wigner crystal in the real space. 



IV. HIGH CONCENTRATION OF 
POLYELECTROLYTE. WEAK SCREENING 
CASE. 



In this section, we deal with the case when Up > Us/Nq 
and there is shortage of spheres, each PE cannot get the 
optimal number, A^o, of spheres found in previous sec- 
tion. In this case, the number of spheres per PE is fixed: 
N = Us /rip. Therefore 



Usq/rip 



Tlp' 

n„ 



S 
C 



(19) 



SO that PECI becomes weaker and linearly decreases with 
Up — Upi as Up grows. When Up increases beyond the 
isoelectric point ripi = UsQ/Q, the total charge of the 
complex Q* is negative. The ratio Q*/Q increases lin- 
early from zero and eventually saturates at unity as Up 
increases further. The behavior of Q^/Q as function of 
Up is plotted by the solid curve in Fig. 2. 

Let us now discuss the behavior of the net charge of 
the sphere = r]{S + x) as Up increases. To do so, one 
has to solve Eq. (M) and find the distance x by which 
the spheres are separated along the PE (we stress again 
that we are interested in the complex far enough from the 
isoelectric point, so that x^ R and Eq. (10) is valid.) 

As Up increases bey ond Us/Nq^ the last two logarith- 
mic terms in Eq. (10) are still negligible compare to the 
second term. Therefore, x is given by Eq. (hsh (it should 
be noted that, here, 6 = C{l — np/npi) is a given length). 
Correspondingly, decreases linearly with 5. 

As Up moves closer to the isoelectric point, the net 
charge g'* = iri{5 + x) decreases. When 



5 < 5 c 



Rln 



L/N - 



R 



R\n{C/R) 



(20) 



(here, we replace ln((L/A^ — x)/R) by \n.{C/R) because 
near the isoelectric point, 5,x ^ C ^ L/N) the fourth 



and the first terms of Eq. (|10|) start to dominate over 
the second and third terms. This gives 



X c^i \S\^ 



\n{C/R) 



(21) 



Of course, Eqs. (H) and (|l|) match each other at 
S = Sc. To continue, let us consider the two important 
limiting cases InA/" <C ln{C/R) and InA/" > ln{C/R). 



A. The case In A/" < \n{C/R) 

In this case, x <C \S\ and, therefore, the charge of a 
sphere r]{S -\- x), decreases to zero and becomes negative 
as Up passes through Upi (see Fig. 2). At Up > Upi^ this 
SCI is driven by the fourth term of Eq. ([l^) : the correla- 
tion energy of PE segment at the surface of the spheres. 
The charge inversion ratio S = \S x\/C |1 — np/npi\ 
(see Fig. 2). 

As Up increases further, the charge of the sphere S -\- x 
grows and the second and the fourth terms of Eq. ( p!o| ) 
become the dominant ones. This gives 



(7*/^ = + x = -Rln 



L/N - 



R 



-i^ln-, (22) 



so that the charge inversion ratio reaches its maximal 
possible value (see Fig. 2) which is equal to that for jthe 
complexation of a single sphere and a polyelectrolyteE3: 



S = 



q 



R^ C 
^ In — . 
C R 



Therefore, roughly speaking, it is inversely proportional 
to the number of turns of PE around the sphere. 

As Up continues to increase, x increases and the third 
term of Eq. ( p^Oj ) becomes important making the sphere 
charge less negative. When x > the second and third 
terms of Eq. ( p!Q| ) dominate. This gives 

L 



mn—, (23) 



SO that the the net charge of the sphere changes sign from 
negative back to positive (not shown in Fig. 2). However, 
the condition of low salt solution, x < vd^ assumed in the 
derivation of Eq. (^oj) , makes this re-entrant inversion of 
charge unrealistic. In practical situation, < C so that 
the necklace remains in the SCI range. A detail consid- 
eration of the strong screening case ro < x is presented 
in the next section. 



B. The case In A/" > \n{C/R) 



In this case, Eq. ( |2l| ) gives x \S\ and the charge 
of a sphere r]{S -\- x) touches zero but stays positive as 
Up passes through the isoelectric point despite the fact 
that the total charge of the complex Q* = N5 changes 
from positive to negative. This is because for a long PE, 
the macroscopic energy is very large and the complex is 
under a strong stress to increase x in order to reduce 
this macroscopic energy. This decreases the amount of 
PE that can wind around each sphere making the sphere 
positive. 

As Up increases beyond the isoelectric point, \S\^ x and 
the net charge q* = S -\- x of each sphere increase as 
well. Eventually q* c:^ q, the PE unwinds from all of its 
spheres and becomes a straight rod to which N spheres 
are attached to. Substituting x = L/N into Eq. ( p!o| ) 
and neglecting the last two terms of this equation, one 
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can estimate the value L/N ai which PE unwinds from 
the spheres: L/N 



mcreases 



C{1 + ^C/R\nM). As 
further, /q saturates at unity. The behavior of as 
the function of Up is depicted by the dashed line in Fig. 
4. 



-1 



q'^/q 




Q*/Q 




/ r 


so / 


V 



FIG. 4. Schematic plot of / q and Q* /Q as functions of 
PE concentration at a fixed and large sphere concentration 
Us for the case B, Sec. IV. There is no SCI in this case. 
The shaded stripe shows the range of Up around Upi where 
condensation of PE molecules happens. 

We would like to emphasize that the inequality 
RhiM ^ \n{C/R) may require unreasonably large 
screening radius ro, so that the behavior presented in 
Fig. 2 for case A of this section is more generic. 



V. CONDENSATION OF PE-SPHERES 
COMPLEXES NEAR THE ISOELECTRIC POINT. 



Now, let us discuss properties of the system near 
the isoelectric point. Exactly at the isoelectric point 
Up = ripi^ the spheres-PE complex is neutral, Q* = q* = 
6 = X = {) and L/N = C. From Eq. one gets the 
energy of one complex as Lhi{A/a). It is the self energy 
of the PE Lhi{R/a) (the PE is straight up to distance 
R) plus the correlation energy —Lhi{R/ A) gained by ar- 
ranging PE turns into one-dimensional Wigner cr yst al at 
the sphere surface (see the discussion after Eq. (p^). A 
consequence of this interpretation for the energy is that 
at the isoelectric point PE molecules condense onto each 
other forming a macroscopic neutral bundle. This is be- 
cause the density of PE in the region where the spheres 
touch each other (the region bounded by broken line in 
Fig. 5) is doubled. Thus, the distance between PE seg- 
ments, At^ is halved. At = A/2, which results in a gain 
in the correlation energy. Simple geometrical calculation 
shows that this region has the area AR. Therefore, the 
PE in this region has total length R. The correlation 
energy gain per unit volume is 




FIG. 5. Cross section through the centers of two touch- 
ing spheres with worm-like (gray) PE wound around them. 
At the place where two spheres touch each other (the region 
bounded by the broken line) the density of PE and the back- 
ground surface charge doubles which in turn leads to a gain 
in the correlation energy of PE segments. Near the isoelec- 
tric point, this gain is responsible for the condensation of 
spheres-PE complexes, forming a large neutral bundle of PE- 
spheres complexes. 

Because of this finite gain in the correlation energy, 
there is a finite range of Up around Upi that the PE 
molecules are still in a condensed state. Let us try to 
find the width of this region. 

To find the boundary of the condensation region on the 
left side of the isoelectric point (up < Upi), one needs to 
compare the total energy of the system in the condensed 
and dissolved states. Here, the condensed state contains 
a macroscopic neutral bundle of PE-spheres complexes 
and {us — UpL/C) leftover spheres per unit volume. The 
bundle is neutral because charging a macroscopic body 
costs a lot of energy. The dissolved state is a solution of 
Up isolated PE-spheres complexes per unit volume, each 
PE adsorbing Us/rip spheres. At the condensation con- 
centration Up = Upi^ we have to balance the correlation 
energy gain AEcorr with the loss in the self energy of 
{us — ripL/C) left-over spheres when they change from 
almost-neutralized spheres at the PE molecules to bare 
spheres in solution. Therefore the equation for the con- 
densation point Up = Upi is 



jC 



R= n 



2R 



or 



^ Upi ^ R^ 



(24) 



(25) 



On the right side of the isoelectric point [rip > Upi), the 
condensed state is a macroscopic neutral bundle of PE 
and spheres and {up — Upi) leftover bare PE molecules. 
AEcorr needs to be balanced with the lost in the self- 
energy of PE molecules when they change from almost- 
neutralized state to bare state in solution. This gives for 
the resolubilization concentration n^^: 



"-^-i^lnA 



At 



ripL^ 



(ri; 



pr 



npi)Lln — 



(26) 
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or 



R 



£ln(ri)/a) 



Finally, the total width of the region, Arip = rip 
around n^i where condensation occurs is 



(27) 



£2 



R 



R 



C\n{rD/(i) £ln(ri)/a) 



Comparing this with Eq. (|2Q|), we see that the width 
of the condensation region is small, well within the re- 
gion where the correlation energy (the fourth term in 
Eq. (|Tq|)) is important in determining conformation of 
the system. Therefore, in this range AEcorr indeed dom- 
inates all other energies in Eq. (|^) as we assumed. 



VI. STRONG SCREENING BY MONOVALENT 
SALT. 

Until now, we assumed the salt concentration is small 
enough so that the screening radius is larger than the 
distance between neighbouring spheres, x. In the case 
of higher salt concentration when td <^ our theory 
needs some modifications. First, the macroscopic energy 
term (the first term in Eq. (^)) has to be replaced by the 
sum of repulsion energies of neighbouring spheres. When 

<C r/) <C X, it still has the form of interaction of two 
point-like charges: 



F{N,x) =N 



q 



* 2 



x^2R 



-{x+2R)/rD 



^Nf{x) . (29) 



At the same time, all the logarithmic factors in Eq. (|^) 
for fix) are cut off at vd instead of x. Correspond- 
ingly, Eq. ( P!q|) (which is the result of the minimization 
of F{N,x) with respect to x at a given N) should be 
replaced by 



dF 

dx 



{S + xf 



-x/ro 



XVd 



R 



5 + x 
R 



. 



In 



R 



(30) 



Let us concentrate on the FECI regime when Up < 
Tig/Nn. In this case, the last two logarithmic terms in 
Eq. (|3Q|) can be neglected. This gives 



td In 



{S^x)R 



, CR 
vd In^- 



Thus, the condition td <C x is equivalent to r/^ <C yCR. 
In this case ri)<Cx<C^~>C. One can see that x only 
weakly depends on the number N of spheres attached 
to the FE. This is because the macroscopic self-energy 
of the complex which forces the FE to unwind from the 



sphere is strongly screened and diminished at distances 
beyond <C x. 

Substituting Eq. ( ^ back into Eq. ( p9|) and opti- 
mizing F{N^x{N)) with respect to we have, to the 
leading term 



{L/N)^ _ rD((5 + X + 2L/N) Cx 



2R 



R 



R 



(28) so that 




and the charge inversion ratio is 



V = 



N6 



C/rn 



HCR/rj,) 



(32) 



(33) 



(34) 



Comparing these results with those of Sec. Ill we see 
that due to screening, the spheres are closer to each other 
{x (X td instead of C}I^R}I^) and a smaller length of the 
FE is wound around a sphere. In other words, the posi- 
tive net charge of each sphere is larger [Lj^ cx C^l^ry^^ 
instead of C?I^R^I^). Therefore, more spheres are at- 
tached to the FE, making charge inversion much stronger 
(J^ (X C}l^ instead of C}l^). At the same time, when 
tb increases to about \J CR^ x ^ L/N C?I^Ri^l^^ 
V ^ (^CjK)^!^ and we come back to the weak screening 
case. 

It should be stressed that, for optimization with re- 
spect to X, the gain in a sphere's self-energy when FE 
winds around the sphere (the second term on the right 
hand side of Eq. (|3Q|)) is balanced with the repulsion 
from its neighbouring spheres (the first term). However, 
when determining and V from Eq. (|3^), the repulsion 
between the spheres described by the first term on the 
right hand side, which is of the order Cvd/R^ is negligible 
compared to the second term Cx/R. This term originates 
from the fact that when one brings a sphere from solu- 
tion to the FE, hence gains the self-energy {L/N)'^ /2R^ 
the FE unwinds from other spheres in order to prepare 
the linker x for this new sphere. 

Let us now consider even smaller screening radius 
vd ^ R' In this case, one has to modify all the energy 
terms of Eq. (29). The self-energy of each sphere be- 



comes (f^TBl2W instead of q"^ /2R and the interaction 
between neighbouring spheres is 



{q*TllR^f 



-x/ro 



(31) As a result, the minimization with respect to x gives 



Now, an equation similar to Eq. (p3) gives 



(35) 



(36) 
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and 




(37) 



SO that the charge inversion is indeed str ong er in this 
case and increases even faster than Eq. (34) with de- 
creasing Equations (|33) and ( |37| ) match each other 
when rjj ^ R. Equation is their combination. 

When the screening length becomes smaller than 
R^/C^ the logarithmic factor ln{CrD/R^) should be re- 
placed by unity and Eqs. ( [35| ) and ( p6| ) give L/N ^ R 
and X <^ R. This means that PE is a straight rod with 
the bare spheres closely packed on it. The number of 
spheres attached to the PE reaches its maximal possible 
value = A^max = L/R^ and so does the charge inver- 
sion ratio V = Vmax = C/R. The behavior of 7^ as a 
function of the screening length td shown by the solid 
line in Fig. o. 
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FIG. 6. Schematic plot of the charge inversion ratios V (the 
solid line) and Vstm (the dashed line) as function of screen- 
ing length ro. Vstm > V td <. C?'^R^'^. V saturates 
when TD R? I a while V^tm saturates when td ^ R- For 
the definition of Testis, see Sec. VII. 

It should be noted that at very small value of the 
screening length, when the energy of interaction between 
a sphere and the PE is less than /c^T, the spheres detach 
from the PE and V rapidly decreases to zero. 

Until now we have concentrated on the effect of strong 
screening on the optimal configuration of the PE-spheres 
complex at small concentration rip < Us/Nq^ when 
spheres are in abundance. Now we want to discuss the 
role of screening in the opposite case, Up > Us/Nq. In 
this case, the number of spheres per PE N = ris/rip is 
fixed and Eq. ( p!9| ) remains valid (screening does not af- 
fect Q*, because in any case all spheres are adsorbed by 
PE.) Qualitatively, Fig. 2 remains valid in this case. The 
length X is now given by the first equality of Eq. jsi] ) for 
td > R and by Eq. ( |35| ) for < (the second equal- 
ities in these equations is not valid because ^ is a given 
length). Therefore, x decreases as S decreases. 

In the case ro > R, when S decreases below r'j^/R, 



Eq. ( pH ) gives x < td and we come back to the weak 
screening case described in Sec. IV, case A. All discus- 
sion about SCI and condensation of complexes in Sec. 
IV, case A remains valid in this case. 

On the other hand, in the case A < vd < R^ our the- 
ory needs some correction. The value of 5 below which 
the correlation energy between PE turns at the surface 
of a sphere is important is given by 



VD R 



\n{C/R) 



(38) 



instead of Eq. ([2Q[). This is because, in this case, the sec- 
ond term in Eq. (|1C)| ) is {S-\-x)rD / R^ instead of {S-\-x) / R. 

At the same time, SCI effect at Up > Upi is strongly en- 
hanced in a way similar to the case of one spherdlJ. This 
is because the charging energy cost for SCI is strongly 
suppressed at small ro while the short-range correlation 
energy between PE turns responsible for SCI remains un- 
affected. At small screening length, S can be larger than 
unity. 

Strong screening {A <^rD <^ R) also affects the range 
of Up where PE molecules form neutral macroscopic bun- 
dle. When Td <C i?, the self-energy cost jC'^/2R m 
the right hand side of Eq. (|2^ ) has to be replaced by 
C'^rn/R^ per sphere while the short-range correlation 
energy on the left hand side remains unaffected. This 
increases the width of this region to 



R^ 



R 



C^td C\ii{rD/o) 



(39) 



This width continue to grow with decreasing vd- When 
rn ^ R^/C - A, R^/C^vd ^ R/C and the width more 
than doubles. When < A^ one can neglect the sec- 
ond term of Eq. ( |39| ) and arrive at Eq. (^. This equa- 
tion predicts a strong growth of Aup/upi with decreasing 
ri), in qualitative agreement with experimental results of 
Ref. |l|. It should be noted again that, as one see from 
comparison with Eq. ( |3^ ) , this coaservation range is well 
within the region of S where the correlation energy be- 
tween PE turns is the dominant energy term. 

Finally, when Up ^ n^^, there is very small number 
of spheres per PE such that the length of the PE linker 
between them is larger than the optimal x given by Eq. 
(|l|) and (H), the linker is no longer straight and each 
sphere with PE wound around it behaves independent 
from each other. SCI saturates at that. given for the case 
of one sphere - one PE complexatior£3. 



VII. POLYELECTROLYTE WITH EXTREMELY 
LARGE PERSISTENCE LENGTH 

In this section we assume that the PE has an extremely 
large persistence length such that it cannot wind around 
a sphere. In this case, the PE is a straight rod to which 
the spheres are attached to. We are interested in the 
PECI regime where the concentration of spheres is large. 
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For a rod-like PE, x = L/N, Nx' = -L/N = -x, 
S -\- X = C In the case of weak screening, vd ^ L/N^ 
Eq. (|l^) can be rewritten as 



r—lnATc^ rin^-— 



(40) 



The physical meaning of this equation is very simple: the 
left hand side is the macroscopic charging energy cost 
when a sphere is brought from the bulk solution to the 
PE. The right hand side is the gain in the correlation en- 
ergy of the Wigner crystal of spheres along the PE which 
helps to overcome the charging energy cost. This corre- 
lation energy is the interaction of the sphere with two PE 
tails of length L/2N which forms a Wigner-Seitz cell. 
The charge inversion ratio can be easily calculated 



N6 _ \n{L/2NR) \n{C/R) 



In A/" 



In A/" 



(41) 



In the case of strong screening, vd <C L/N ^ the macro- 
scopic charging energy cost should be replaced by the 
repulsion between neighbouring spheres. At the same 
time, the logarithmic term in the expression for the cor- 
relation energy of the Wigner crystal of spheres along the 
PE should be cut off at td instead of L/N. Eq. ( ^ ) now 
reads: 



£2 

— ( 



-^/^'■^ =£111 



R 



(42) 



which gives L/N :^ ri:)ln(£/ri:)), and the charge inver- 
sion ratio 



c 



vd \n{C/rD) 



> 1 



(43) 



Let us now compare these results with those for an 
intrinsically flexible PE case studied in Sec. Ill (weak 
screening) and Sec. V (strong screening). 

At weak screening (Sec. Ill) vd > C > R\nJV^ Eq. 
(|l^ and Eq. ^ give 



T^stiff 



R\n{C/R) 



\n{jC/R) 



£1/4(7^ In Ar)3/4 (£/i^)i/4(inAA)3/4 



< 1 



The last inequality is due to C/R^ 1 and hiM ^ 1. 

vd decreases further so that C > vd > VCR we 
enter the strong screening regime for the rod-like PE but 
still stay in the weak screening regime for flexible PE. 
Using Eq. ( p3| ) for T^stiff and Eq. (|l^) for V , we can 
easily see that the charge inversion for the rod-like PE 
starts to become stronger than charge inversion for the 
flexible PE when vd ^ C?I^R^I^. 

Wh en td continues to decrease in the range td < 
\/ CR^ we are in the strong screening regime for both 
types of PE. Equation ( |3^ ) and ( |37| ) show that V is 
of the order of \/ C/vd which is much smaller than 
Vsw ^ C/td- 



The behavior of the charge inversion ratio as a func- 
tion oi td for the two types of PE is shown in Fig. ||. 
One can explore a transition from the flexible PE case 
to the stiff PEp^se by applying an external stretching 
force to the PEEII. To describe this phenomenon, one has 
to add to the free energy (|^) of the complex an addi- 
tional term —J^N{x-\-2R), where is the external force. 
This new term is linearly proportional to the length of 
the spheres-PE complex. It adds a negative term —J^ to 
the right hand side of Eq. (p!o| ) and therefore increases x. 
One then can proceed in exactly the same way as in Sec. 
Ill to find the conformation of the complex. At weak 
screening when vd ^ C?I^R}I^ (Sec. Ill), it is not diffi- 
cult to show that x increases linearly with the strength of 
the external force when this force is small. At the same 
time, A^o decreases linearly with T so that one by one, 
the spheres leave the PE as the force increases and PECI 
becomes weaker. When T ~ {rf / D){C/ R — In AT) (so 
that the force helps to balance the attractive potential 
of the sphere with the macroscopic repulsive potential of 
the complex), one obtains x ^ L/N and the PE unwinds 
completely from the spheres and becomes straight. The 
problem of complexation of PE and spheres, then, be- 
comes that of a stiff PE described at the beginning of 
the section. This sequential release of spheres is simi- 
lar to thjC problem of stretching a PE necklace in poor 
solventE3. 

The picture is completely reversed at the strong screen- 
ing case when td <C C?I^R}I^. In this case, as one 
stretches the PE, the spheres come to the PE one by one 
and make the charge inversion stronger. The strength 
of the force at which PE unwinds completely from the 
spheres can be calculated in exactly the same way as in 
the weak screening case. In the strong screening case, 
however, the macroscopic repulsive potential of the com- 
plex is very small so that the external force has to over- 
come only the attractive potential of the sphere in or- 
der to unwind the PE molecule. Therefore, PE un- 
winds completely when T ~ rfC/RD for > R and 
J" - r]^CrD/R^D for vd < R. 

It would be interesting to verify experimentally that 
the spheres leave the PE at weak screening and condense 
on the PE at strong screening when the PE undergoes 
an external stretching force. 



VIII. CONCLUSION 

In conclusion, we would like to discuss four most im- 
portant approximations used in this papers. Let us start 
from the use of the Debye-Hiickel linear theory to de- 
scribe screening by monovalent salt. It is known that if a 
PE molecule or a sphere are charged strongly enough this 
linear approximation does not work and the nonlinear 
condensation of counterions takes place, which leads to a 
renormalization of their charge. For the case of a rod-like 
PE this phenomenon is known as the Onsager-Manning 
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condensationE2l. It happens when the hnear charge den- 
sity of PE T] is larger than r]c = DkBT/e^ where /cb is 
the Boltzmann constant and T is the temperature. Cor- 
respondingly, Debye-Hiickel theory used above is valid 
when T] < ric 

The condition for the absence of counterion conden- 
sation on the charged spheres is more complicated and 
involves the concentration of monovalent salt as well. It 
is known that if a sphere is charged strongly enough, its 
counterions condense onto its surface to reduce its charge 
to the universal critical value qc = DRekBT\n{cs/c)^ 
where e is the elementary charge, Cs ~ qc/R^e and c 
are the counterion concentrjations at the sphere surface 
and in the bulk respectivelyEll. The condensation on the 
spheres can be neglected \i q = Cr] is less than qc or 

ro > i^e^^/^^^ . (44) 

At small enough r] (such that r]/r]c < R/^)^ this condi- 
tion reduces to > R. 

When rjj < R^ each sphere can be considered as a 
plane with surface charge density q/AiiR^ and it is also 
known that li is small enough, a charged plane is lin- 
early screened. Specifically, Eq. (73) of Ref. ^ shows 
that screening is linear if 

vd < Ae'^^l'^ - 7"^"^'^"^' ^^^^ 

At small enough r] (such that r] < rid ln(>C/i?)), this con- 
dition reduces io vd < R- 

Thus our theory has a wide range of applicability. For 
a PE with 

T] < rjcR/C (46) 

it is applicable for any value of the screening length 
Remarkably, the same inequality for r] also guar- 
antees that no Onsager-Manning condensation occurs on 
the spheres-PE rod-like complex with the inverted linear 
charge density r]V^ even though the magnitude of this 
inverted charge can be much larger than the bare charge 
of the PE {V ^ 1). Indeed, as we already know, C/R 
is the maximal possible value for V (see Fig. 6) there- 
fore Eq. ( ^ guarantees that rjV is smaller than the 
Onsager-Manning critical linear charge density r]c. 

On the other hand, our theory literally is not applica- 
ble to strongly charged PE such as double helix DNA or 
too strongly charged spheres. We do, however, believe 
that our main results are qualitatively applicable in this 
case with properly normalized charges of the spheres and 
PE. A more detail study of the condensation problem will 
be the subject of a future work. 

The second approximation which we want to address 
here is the assumption that PE is fiexible. It is valid 
if the elastic energy of PE winding around a sphere is 
smaller than the Coulomb energy of complexation. This 
can be true even if the persistence length of the strongly 
screened PE is of the order of a sphere diameter or some- 
what larger. For example, theoretical estimates show 



that even in the nucleosome bead, the rigidity of DNA 
plays a secondary role. We should emphasize here that 
we talk about the intrinsic persistence length, ignoring 
the well known Odijk-Skolnick-Fixman electrostatic en- 
hancement of the persistence length, which can be very 
substantial for a free PE pplecule in the solution with 
small concentration of saltEj. The reason for this is that 
at the sphere surface the charge of PE is screened by 
the sphere's positive charge and it is the intrinsic per- 
sistence length of PE which determines whether PE can 
wrap around a sphere. All electrostatic interactions are 
explicitly taken into account in the correlation picture 
considered. 

If the PE intrinsic persistence length is much larger 
than i?, the ground state of the complex can strongly dif- 
fer from that in the fiexible case. In Sec. VI we studied 
the extreme case when the persistence length of the PE 
is infinite and the PE is rod-like. There is a wide range 
of intermediate magnitudes of persistence length which 
is not studied here. In this range nontrivial star-like con- 
figurations become possible even for one sphereE3. For 
many spheres one can imagine different kinds of struc- 
tures. Near the isoelectric point, it can be a rod-like com- 
plex made of a PE solenoid densely stuffed by spheres. 
It can be a similar cylinder where PE is instead winding 
around spheres makes several parallel to the cylinder axis 
straight strands on the surface of the cylinder, which are 
connected by loops at the cylinder edges. This strands re- 
pel each other and form periodic in polar angle "Wigner 
crystal" . These and other possible configurations should 
be studied in future works. 

The third simplification used in this paper was an as- 
sumption that the number and the size of spheres in the 
solution is fixed. When we try to compare this theory 
with the data on the micelles-PE systemtl we immedi- 
ately see that, strictly speaking, this is a different prob- 
lem, because in the experiment of Ref. [l] the amount of 
lipids is fixed, but the number and size of micelles is de- 
termined by equilibrium conditions and may depend, for 
example, on the screening radius of the solution. This 
complication should be taken into account by a future 
theory. 

The fourth important assumption we made in this pa- 
per is that the concentration of spheres, in the solution, 
ng, is large enough so that their entropy can be neglected. 
This assumption leads, for example, to the conclusion 
that in the case Up > Us/Nq, all the spheres are consumed 
by PE. Actually, even in this case, there is a finite, but 
exponentially small concentration, n^o, of free spheres in 
the solution next to PE because the binding energy of a 
sphere to PE is finite. In other words, this is the concen- 
tration of the "saturated vapour" of spheres right above 
the correlated liquid of spheres on PE. When total con- 
centration of spheres, n^, is so small that it is comparable 
to Uso effect of free spheres becomes very important. For 
example, at such small rig spheres fail to neutralize PE 
near isoelectric point. At the limit of ^ neutral- 
ization happens at Us = riso- Therefore, the plot of the 
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line ns{np) at which neutrahzation takes place deviates 
from isoelectric line ris = ripL/jC at small Up and (see 
dashed line at Fig. 7). Correspondingly, the region of 
the condensation of PE-spheres complexes on the plane 
(77,5, Up) looks as the shaded region shown in Fig. 7. 
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FIG. 7. Schematic phase diagram of the condensation of 
spheres-PE complexes in the (n^, n5)-plane for a very large 
PE length L. Complexes condense in the shaded domain. The 
straight line corresponds to the isoelectric point ris = UpL/C. 
The dash line shows the function ris (rip) at which an isolated 
PE-spheres complex is neutral. 

Theory of condensation at small rip is similar to the- 
ory of condensation of DNA biL-multivalent counterions 
at small concentration of DNAllZl. Experimental data on 
re-entrant condensation of DNA have-qualitatively simi- 
lar shape of the condensation domainEO. 

Summarizing our results, we have studied complexa- 
tion of a negative PE with positive spherical macroions 
in salty water. Under solely the influence of electrostatic 
interactions, a PE molecule winding around individual 
spheres binds spheres in a beads-on-a-string structure. 
At a large sphere concentration, we found interesting 
phenomena on both sides of isoelectric point, at which 
the total charge of all PE molecules is exactly compen- 
sated by total charge of spheres. When the PE concen- 
tration is below the isoelectric point, spheres overcharge 
PE so that the net charge of PE together with bound 
spheres is positive and can be substantially larger than 
the absolute value of the bare charge of PE. When the 
PE concentration is above the isoelectric point, the net 
charge of PE is not inverted, while the net charge of each 
sphere together with the PE segment winding around it 
becomes inverted and negative. It can be larger than the 
bare charge of the sphere if the screening radius of the 
solution is small enough. 

In the narrow vicinity of the isoelectric point PE- 
spheres complexes condense together in bundles. We cal- 
culated the width of the range of the PE concentration 
where the condensation takes place and showed that it 
grows very fast with decreasing screening radius of the 
solution. 

All these phenomena are results of repulsive correla- 
tions between PE turns on the surface of spheres and 
of the spheres on PE, which can not be included in any 
Poisson-Boltzmann-like description of the mutual screen- 



ing of PE and spheres. The repulsive correlations be- 
tween PE turns on spheres are responsible for charge 
inversion of individual spheres (SCI) and the condensa- 
tion of PE-spheres complexes in bundles. In the latter 
case this phenomenon is illustrated by Fig. 5. In the 
former case, additional discussion of the physics of this 
phenomenon can be found in Ref. ^ Wigner-crystal-like 
correlations of turns mean that each turn is surrounded 
by a stripe of positively charge macroion surface, which 
can be considered as its positive correlation hole. 

The most interesting phenomenon of the inversion of 
charge of PE (PECI) by a large number of adsorbed 
spheres is related to the fact that each sphere is bound to 
several turns of negatively charged PE. These turns can 
be considered as an correlation hole, because this is the 
part of PE, which almost exclusively interacts with the 
given sphere. In this sense, once more we are dealing with 
correlations. The segments of PE wound around each 
sphere have the same length and they constitute most of 
the PE's total length. Therefore, these correlations are 
similar to Wigner-crystal-like correlations which-aDe-|p.- 
sponsible for charge inversion of a rigid macroionlij t2lE3. 
This confirms our point of view that Wigner-crystal-like 
correlations are the universal driving force of charge in- 
version. 

Our theory, with minor modifications, can also 
describe the complexation of polyelectrolytes with 
macroions of non-spherical shape. An example of 
such system can be the complexation of two oppositely 
charged polyelectrolytes in water solution. Let us assume 
that the negative PE is long and flexible while the posi- 
tive PE is shorter, stronger charged and rigid. Then the 
negative PE molecule wraps around the positive one and 
the only change which should be made in our theory is to 
replace the expression for self-energy of charged sphere by 
the corresponding expression for charged rod. The case 
when the positive PE is flexible on the first glance seems 
to be more complicated. However, away from isoelec- 
tric point positive PE is overscreened or underscreened 
by the wrapping negative one so that it has a rod-like 
shape. Therefore, even in this case our theory is valid 
with the above-mentioned minor change. 
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